van der Waals interaction in nanotube bundles: consequences on vibrational modes 
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We have developed a pair-potential approach for the evaluation of van der Waals interaction between 
carbon nanotubes in bundles. Starting from a continuum model, we show that the intertube modes 
range from 5cm~^ to 60c7n~^. Using a non- orthogonal tight-binding approximation for describing 
the covalent intra-tube bonding in addition, we confirme a slight chiral dependance of the breathing 
mode frequency and we found that this breathing mode frequency increase by ~ 10 % if the nanotube 
lie inside a bundle as compared to the isolated tube. 
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Research in carbon nanotubes is now a very active field 
both because of their fascinating cylindrical structure 
and potential applications (see ||l|,|[). On the characteri- 
zation side, the direct measurement of atomic structure is 
now possible by means of scanning tunneling microscopy 
(STM) experiments P|. This technique gives informa- 
tion on both the local atomic structure as well as elec- 
tronic properties. However it is too time-consuming for 
routine sample characterization. Thus one has to resort 
to more macroscopic techniques, which are able to han- 
dle the whole as-grown sample. One of the most pop- 
ular alternative experimental tools is the measurement 
of the vibrational spectrum by Raman spectroscopy (see 
e.g. Q). Since the suggestion was made Q that the vi- 
brational frequency of the fully symmetric Ag breathing 
mode (BM) of an isolated nanotube could be used to 
determine the nanotube diameter, given its marked fre- 
quency dependence on the diameter of the tube, this idea 
has been widely used as a characterization tool. But, it 
should also be noted that a small chiral dependance of 
BM frequency was recently found by ab-initio calcula- 
tions |^,||. 

However for developing useful technological applica- 
tions, it is of primary importance to refine the character- 
ization techniques of nanotube samples, taking into ac- 
count the local environment in which the nanotubes find 
themselves. Surprisingly, theoretical studies of the reso- 
nance frequency have largely focused on isolated tubes, 
and no such studies have been reported for bundles of 
tubes, which is the form in which the majority of the 
nanotubes are to be found in experimental samples. The 
aim of the present communication is to show how the 
tube-tube interactions can be taken into account for the 
characterization of the BM. Due to its radial character, 
the BM is likely to be the most influenced mode by the 
nanotube packing and therefore it is neccesary to quan- 
tify this effect. But before we address this question, we 
will present an evaluation of the inter-tube vibrational 
modes, which could also allow the characterization of 
samples by inelastic neutron scattering experiments. 



The theoretical model we have developed to deal with 
the tube-packing effects in the vibrational spectra com- 
bines a description of the covalent bonds inside each nan- 
otube via a non-orthogonal tight-binding parametriza- 
tion 0], which has proved to work very well for the 
structural and mechanical properties of carbon nan- 
otubes [p|, with a pair-potential approach to deal with 
the van der Waals interaction between carbon nanotubes. 
The long-range dispersion interaction is described by 
a carbon-carbon Lennard- Jones potential Vcc given by 
Vcc — ~^ + §rT J where d is the carbon-carbon distance 
and Cg and C'12 are constants fitted to reproduce the 
structural properties of graphite [gj . Despite its intrinsic 
simplicity, the present approach is widely used to simu- 
late van der Waals interactions between molecules (see 
PQI ) . In all the calculations described below we consider 
infinitely long tubes. Let us also note that, although 
ab-initio descriptions are in general more accurate and 
reliable, they are very computationally demanding and 
are limited in the size of the system PJll| . 

First, to study the intertube modes, we simplify further 
the model and we consider carbon nanotubes as continu- 
ous cylindrical surfaces of density a. In this case we can 
compute analytically the potential felt by a carbon atom 
situated at a distance a from the center of the nanotube 
of radii R. The potential reads {a > R): 



VMa) - 7vri?a ( ^F [5/2, 5/2, 1, (R/a)' 
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+ |§^^[llAll/2,l,(i?/a)2]| , (1) 

where F is the hypergeometric function [|l2[ . This model 
cannot distinguish between different layer-stackings in 
multilayer compounds, nor can it take account of the 
tube chiralities when applied to bundles of nanotubes. 
We feel that this is a good approximation as the registry 
property in nanotube bundles is frustrated by the geom- 
etry of the bundle. Our model thus can be considered 
to represent an average over the different stacking pos- 



sibilities in a bundle of tubes of similar diameter havir 
chiralities compatible with that diameter. In this sens 
we expect the continuous model potential to be able 1 
give a realistic description of the intertube modes. Th 
approach is similar to the one developed by Girifalco ft 
the Ceo fuUerenc |13|. Within this simplified model tl 
van der Waals interaction energy per unit length betwee 
two carbon nanotubes, Vtt, follows immediatly by ni 
merical integration of Vcria) over the surface of the sei 
ond nanotube. Now, the total energy of a bundle wit 
N tubes, Vb, is obtained by summing Vtt over all tubi 
pairs. Before the computation of the vibrational spectr; 
the total energy for a given bundle is minimized wit 
respect to the tube positions JlJ]. Then, the intertul 
vibrational modes (eigenvalues and eigenvectors) are di 
duced from the diagonalisation of the dynamical matr: 
constructed from a numerical evaluation of the secoE 
derivatives of the total energy with respect to the tut 
center coordinates |15|| . As a practical remark, we ha'\ 
built the finite-size bundles enforcing the 6-fold symm 
try in the section of the bundle, which has one nanotul 
at its center (the number of tubes is A^ = 6i -t- 1 wit 
i being the number of hexagonal-shells in the bundle 
In this case, all the inter-tube vibrational spectra can I 
sorted according to the Cgt, point group. In particular v 
find 3i Raman active modes and 2i infrared (IR) acti"\ 
modes. 

In Fig. ^ we show the Raman and IR active modi 
of a finite size bundle of (10,10) {R—6.8A) armchair nai 
otubes as a function of the number of tubes N. Intertul 
modes are found between bcm~^ to 60cm~^. To the bei 
of our knowledge, no experimental evidence of the obse 
vation of Raman active modes in this energy range have 
been reported. The origin of these lack of observational 
evidences is probably a low scattering cross section and 
difficulties to deconvoluate the elastic-scattering peak. 
Inelastic neutron scattering experiments should also be 
able to probe such modes even if the difficulty of rigid 
displacement of long, massive cylinders could reduce the 
excitation probability. However such experiments require 
a large quantity of highly purified nanotube powder and 
at present no data is available for the inter-tube vibra- 
tion energy range llq] . For comparison between our sim- 
ulation and future experimental work, we show by the 
solid curves in figure 3 the total vibrational density of 
states of a bundle consisting of 55 tubes. Note also that 
the squashing mode of isolated tubes and the libration 
mode are predicted at about 16cm^^ ||l^ and are 
likely to interfere with intertube modes. This last mode 
is expected to play a role in temperature dependant con- 
ductivity of nanotube bundles since it can induce temper- 
ature dependent hopping conductivity. The intertubes 
mode described here could have similar behaviour. 
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FIG. 1. Intertube mode frequencies of Cgu bundle of in- 
finitely long nanotubes of radius 7?=6.8A as a function of the 
number of tubes. Only Raman active (Ai and E2) modes and 
IR modes (Ei) are represented by fill/open circles . The total 
vibrational density of states, using a gaussian broadening of 
0.8cm~^, for a bundle with 55 tubes is also given. 

After presenting the new data on intertube vibrational 
modes, we look at the influence of packing in the BM. The 
calculations were performed in a frozen phonon approach 
using the hybrid tight-binding plus continuous Lennard- 
Jones potential described above. We use a conjugate- 
gradient relaxation scheme to determine the geometry 
of the isolated tube before making the full-relaxation of 
the tube-bundle. After relaxation, we applied the frozen 
phonon approach to evaluate BM frequencies by comput- 
ing the total energy (intra- and inter- tube) for a 0.1 % 
change of the tube radius. One of the limitation of the 
frozen phonon calculation for the BM in bundles is the 
fact that the pure radial mode is no more an eigenmode 
of the bundle system since packing breaks the symmetry 
of isolated tubes. This will lead, in addition to the fre- 
quency shift, to a degeneracy lift. These effect were study 
for Ceo in |19[- Here we do not take this symmetry low- 
ering in consideration but we believe that frozen phonon 
approach catches the essential feature on the packing ef- 
fect on BM modes and gives a first good evaluation of 
the resonance frequency increase that will be usefuU for 
Raman spectroscopy analysis of nanotube sample. 
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FIG. 2. Breathing mode (^i) frequencies as a function of 
the tube radius. This work : Circles are for armchair tubes, 
squares for zig-zag and diamond for chiral tube. Solid symbols 
are for isolated tubes whitin tight-binding approach and open 
symbols are for tubes in a bundle (see text). + are estimated 
from force field approach [E0[, x |q| and * H come from 
ab-initio calculations. 

In Fig. y and Table | we present the results for the 
BM modes of isolated tubes of different radii and chi- 
ralities as well as the values for the corresponding in- 
finite bundles. As far as isolated tube concerns, we see 
that armchair (n, n) and zigzag (n, 0) tubes do not follow 
exactly the same scaling law with diameter. We repro- 
duce the known l/R scaling of the BM frequency {v), in 
particular the fit of our data to z^ = C/R give the fol- 
lowing values: C =1307cm~^ A for armchair tubes and 
C=1282cm~^ A for zigzag tubes. Chiral tubes BM mode 
lie between armchair and zigzag lines. Recent ab-initio 
results H showed C=1180cm-i A and C=1160cm-i A 
for armchair and zigzag tubes, respectively and force 
field data are reproduced by C=llllcm~^ A pO] or 
C=1147cm~^ A pll irrespective of tube chirality. Con- 
sidering that tight-binding methods usually overestimate 
vibrational frequencies by 5-10 % P3| , p4| , the results pre- 
sented above are not surprizing. We emphasize that our 
tight-binding calculation reproduce well the slight dif- 
ference between armchair and chiral tubes found by ab- 
initio methods ^. The reason for this success is con- 
nected with the non-orthogonality of our tight-binding 
method that correctly describe the re-hybridization of 
carbon-carbon bonds when the nanotube breathe. More- 
over, the low computational cost of the TB approach 
compared to DFT allowed us to perform a more system- 
atic study of chiral and achiral nanotubes (see Table |) . 

We can now turn to the case of nanotube bundles. 
As expected, Fig.0 and Table H show a clear increase in 
the BM frequencies when tubes are packed into bundles. 
The relative increase goes from 5 % for a tube of ra- 
dius i?=3.5A to 15 % for R^Sk. For R^G.Sk (that are 
revelant considering actual produced nanotubes), the fre- 
quency shift is evaluated to be of the order of 10 %. Very 



recently, Venkateswaran et al. |25|, using a similar model 
to study the pressure dependance of nanotube bundles 
Raman modes, found a 8% increase of the BM of (9,9) 
tube when they are packed compare to isolated tube. 
This is consistent with the present result and the differ- 
ence (we found 10% increase for a (9,9) tube) is likely to 
come from a sligth difference in parameters used both for 
the tight binding and for the pair van der Waals poten- 
tial. 

Focussing on consequences on the interpretation of 
Raman spectra, the effect of neighboring tubes on BM 
frequencies can lead to misinterpretation of Raman re- 
sults and to a 10 % error in radius determination. Even 
if tight-binding calculations do not give exact numbers 
for vibrationnal energies, our results are inclined to con- 
clude that the experimental Raman peak at 180cr7i^^ is 
to be associated with tubes of diameter larger than the 
(10, 10) tube diameter. Indeed, simulation for single-wall 
nanotubes predicted a BM for a (10,10) tube between 
163c7TT-~^ (force field [gO[) and 194c7tt-^^ (this work) or 
at nScm^^ for the ab-initio evaluation in [^. If we 
consider those results as 10 % under-evaluation of the 
BM of tube in a bundle, we get frequencies higher than 
180cm^^ (except for the less sophisticated force field 
model) and the key frequency of 180cto~^ is associated 
with tube with diameter larger than 6.8A. This conclu- 
sion is consistent with X-rays p6p7| ] and neutron ES] 
diffraction if the tube-tube distance is taken to be 3.2^1 
(and not the graphite interlayer distance 3.35A) since 
the diameter polydispersity required for diffraction spec- 
tra fit [E8[ cmplies the presence of tube larger than the 
(10,10) tube in bundles. Moreover, electron diffraction 
p9[ on single nanotube rope lead to the conclusion that 
6.8A < R < 7.5A and then to a mean radius larger than 
the the (10, 10) radius. 

As the BM concerns, we have only considered bundles 
made of a infinite number of nanotubes and the study of 
finite bundles will obviously lead to a gradual change of 
the BM vibration frequency going from a frequency close 
to the bulk-value for the central tubes (surrounded by 
6 tubes) to a frequency close to the isolated-tube value 
corresponding to tubes at the bundle-surface. This be- 
havior is expected from the short-range tube-tube inter- 
action (at the scale of tube diameter). Note that STM g] 
and electron diffraction p3| experimental studies of car- 
bon nanotubes concluded that all tube chiralities can be 
found in samples with a very narrow diameter distribu- 
tion. This will not change drastically our conclusions, 
even more, if tube chirality are randomly distributed 
whithin the rope, no registry is possible between adja- 
cent tubes and the continuous model is well justified. 
Furthermore, we have checked that an increase/decrease 
of the tube diameter within the rope leads to a reduc- 
tion/increase of the BM vibrational frequency of neigh- 
bouring tubes. Then, both the surface effect and 'defect' 
(larger or smaller tubes) will then lead to a broadening 



of the experimental Raman spectra as observed experi- 
mentaly. 

TABLE I. Breathing Mode (BM) frequencies for various 
carbon nanotubes. The first column defined the tube, the 
second one is the radius of the relaxed structure. The fre- 
quencies of the isolated tubes and tubes in bundle are given 
in the next two columns. The last column is the relative in- 
crease of the BM frequency when bundles are packed. 

{n,m) R{A) v^soi{cm~'^) Vbundie{cm~'^) shift(%) 



(6,4) 


3.45 


(8,2) 


3.63 


(7,4) 


3.81 


(10,0) 


3.91 


(6,6) 


4.07 


(10,1) 


4.16 


(11,0) 


4.34 


(12,0) 


4.70 


(7,7) 


4.81 


(10,4) 


4.92 


(13,0) 


5.11 


(12,3) 


5.41 


(8,8) 


5.49 


(15,0) 


5.90 


(14,2) 


5.93 


(9,9) 


6.17 


(12,6) 


6.23 


(10, 10) 


6.85 


(16,4) 


7.19 


(11,11) 


7.53 


(20,0) 


7.84 


(12,12) 


8.21 



366 
344 
313 
328 
313 
304 
297 
269 
268 
256 
247 
232 
239 
214 
211 
214 
205 
195 
179 
178 
166 
164 



384 
362 
330 
349 
332 
323 
316 
289 
288 
276 
268 
253 
259 
236 
233 
236 
227 
217 
202 
201 
190 
187 



4.8 
5.3 
5.2 
6.2 
6.0 
6.3 
6.4 
7.4 
7.4 
7.9 
8.3 
9.1 
8.6 
10.1 
10.4 
10.1 
10.7 
11.4 
12.8 
12.9 
14.2 
14.4 



In conclusion, we have presented the first study of the 
inter-tube vibrational modes (both Raman and IR ac- 
tive) in bundles of single walled nanotubes and proposed 
neutron inelastic scattering as a experimental test of the 
validity of our empirical model. We have also showed the 
first computational evidence of the packing influence on 
the BM of nanotubes and drawn conclusions on the inter- 
pretation of the experimental Raman spectra and how to 
extract useful experimental information about the nan- 
otube structure and diameter. 

Acknowledgments: We have benefited from fruitful dis- 
cussions with P.Senet, E.Anglaret, S.Rols, J.L.Sauvajol, 
A.Loiseau and Ph.Lambin. This work was supported 
by the TMR contract NAMITECH(ERBFMRX-CT96- 
0067(DG12-MIHT)), the Belgian Program(PAI/UAP 
4/10) and JCyL(VA28/99). L.H. is supported by the 
Belgian Fund for Scientific research (FNRS). 



[1] G. D. M. Dresselhaus and P. Eklund, Science of 
Fullerenes and Carbon Nanotubes (Academic Press, San 
Diego, 1996). 

[2] M. Terrones, W. Hsu, H. Kroto, and D. Wahon, Topics 
in Current Chemistry 199, 190 (1999). P. Ajayan and T. 
Ebbesen, Rep. Prog. Physics 60, 1025 (1997). 

[3] J. Wildoer, L.C. Venema, A.G. Rinzler, R.E. Smalley and 
C. Dekker, Nature 391, 59 (1998). L. Venema, J.W.G. 
Wildoer, S.J. Tans, J.W. Janssen, L.J. Hinne, T. Tu- 
instra, L.P. Kouwenhoven and C. Dekker, Science 283, 
52 (1999). T. Odom, J. Huang, P. Kim, and C. Lieber, 
Nature 391, 62 (1998). 

[4] P. Eklund, J. Holden, and R. Jishi, Carbon 33, 959 
(1995). A. Rao, E. Richter, S. Bandow, B. Chase, P.C. 
Eklund, K.A. Williams, S. Fang, K.R. Subbaswamy, M. 
Menon, A. Thess, R. Smalley, G. Dresselhaus and M.S. 
Dresselhaus Science 275, 187 (1997). 

[5] J. Kiirti, G. Kresse, and H. Kuzmany, Phys. Rev. B 58, 
R8869 (1998). 

[6] D. Sanchez- Portal, E. Artocho, J.M. Soler, A. Rubio and 
P. Ordejo, Phys. Rev. B 59, 12678 (1999). 

[7] A non-orthogonal tight-binding a + tt scheme was used 
using 11 irreducible points over the one-dimensional 
Brillouin-zone of the isolated tube. The hopping integrals 
and the hamiltonian and overlap matrices are tabulated 
as a function of the internuclear distance on the basis 
of first principle DFT calculations and only one and two 
center contributions to the hamiltonian matrix were re- 
tained (see a detailed description in p3,po|). 

[8] E. Hernandez, C. Goze, P. Bernier, and A. Rubio, Phys. 
Rev. Lett 80, 4502 (1998). 

[9] L.A. Girifalco and R.A. Lad, J. Chem. Phys. 25, 693 
(1956). The parameters of the Lenard- Jones potential for 
carbon are 6*6=20 eV A'' and Ci2=2.48 10* eVA". 
[10] J. Israevlachvili, Intermolecular and surf aces forces (Aca- 
demic Press, San Diego, 1992). 



[11] A. Rubio, D. Sanchez-Portal, E. Artacho, P. Ordejon and 
J.M. Soler Phys. Rev. Lett. 82, 3520 (1999). 

[12] M. Abramowitz and I. Stegun, Handbook of Mathematical 
Functions (Dover Publication, New York, 1968). 

[13] L.A. Girifalco, J. Phys. Chem. 96, 858 (1992). 

[14] For an infinite 2D-packing of (10, 10) tubes with radii 
of _R=6.8 A the relaxed intertube distance is found to be 
3.14 A. Note that the value for the interlayer distance be- 
tween two flat carbon sheets within the continuous model 
is 3.27A and for an infinite layer stacking (graphite) is 
3.23A. This distance correspond to the actual interlayer 
distance in graphite 3.34 A. The '^O.lAdifference is a 
measure of the graphene layer thickness. Then, we have 
confirmed the assumption that curvature reduces the in- 
tertube distance in the bundle as compare to graphite by 
more than 0.1 A. 

[15] We note that for the calculation of phonons in Ceo solid, 
the Lennard- Jones carbon-carbon intermolecular poten- 
tial has been succesfuUy used together with a bond charge 
model [|19|. This last part of the potential was added 
empirically to reproduced the relative orientation of Ceo 
molecules in the solid. In the case on nanotube-bundles 
there is no evidence of such orientation and we do not 
include such bond charge model in our calculation. 

[16] For a review on neutron experiments on nanotubes, see 
S. Rols, E. Anglaret and J.L. Sauvajol Submitted. 

[17] Y. Kwon, S. Saito, and D. Tomanek, Phys. Rev. B 58, 
R13314 (1998). 

[18] This libration mode do not come out the present simula- 
tions because of the continuum character of the tube-tube 
potential we have used. 

[19] J. Yu, L. Bi, R. Kaha, and P. Vashista, Phys. Rev. B 49, 
5008 (1994). 

[20] S. Rols, (Private communication). Calculations were 
done within the model described in ref [M (unpublished) . 

[21] R. Saito,T. Takeya, T. Kimura, G. Dresselhaus and M.S. 
Dresselhaus, Phys. Rev. B 57, 4145 (1998). 

[22] V. Popov, V. V. Doren, and M. Balkanski, Phys. Rev. B 
59, 8355 (1999). 

[23] D. Porezag, Th. Frauenheim, Th. Koler, G. Seifert and 
R. Kaschner, Phys. Rev. B 51, 12947 (1995). 

[24] The same scheme for Ceo gives a BM frequency at 
550cm~^ where experimental value is 497cm~^. 

[25] U. D. Venkateswaran, A.M. Rao, E. Richter, M. Menon, 
A. Rinzler, R.E. Smalley and P.O. Eklund, Phys. Rev. B 
59, 10928 (1999) 

[26] C. Journet, W.K. Maser, P. Bernier, A. Loiseau, M. 
Lamy de la Chapelle, S. Lefrant, P. Deniard, R. Lee and 
J.E. Fischer, Nature 388, 756 (1997). 

[27] A. Thess, R. Lee, P. Nikolaev, H. Dai, J. Robert, C. 
Xu, Y.H. Lee, S.G. Kim, A.G. Rinzler, D.T. Colbert, 
G.E. Scuseria, D. Tomanek, J.E. Fischer, R.E. Smalley, 
Science 273, 483 (1996). 

[28] S. Rols, R. Almairac, L. Henrard, E. Anglaret and J.L. 
Sauvajol, Eur. Phys. J. In Press (1999). 

[29] L. Henrard, A. Loiseau, C. Journet, and P. Bernier, Eur. 
Phys. J. (in press). 

[30] C. Goringe, D. Bowler, and E. Hernandez, Rep. Prog. 
Phys. 60, 1447 (1997). 



